import pandas as pd

from pgmpy.estimators import MaximumLikelihoodEstimator 
from pgmpy.models import BayesianModel
from pgmpy.inference import VariableElimination
import csv

lines=list(csv.reader(open('data7_names.csv','r')))
attributes=lines[0]

#attributes=['age','sex','cp','trestbps','chol','fbs','restecg','thalach','exang','oldpeak','slope','ca','thal','heartdisease']

heartdisease=pd.read_csv('data7_heart.csv',names=attributes)

model=BayesianModel([('age','trestbps'),
                     ('age','fbs'),
                     ('sex','trestbps'),
                     ('sex','trestbps'),
                     ('exang','trestbps'),
                     ('trestbps','heartdisease'),
                     ('fbs','heartdisease'),
                     ('heartdisease','restecg'),
                     ('heartdisease','thalach'),
                     ('heartdisease','chol')
                    ])

model.fit(heartdisease,MaximumLikelihoodEstimator)

heartDisease_infer=VariableElimination(model)

query=heartDisease_infer.query(['heartdisease'],{'age':20})
print(query['heartdisease'])

query=heartDisease_infer.query(['heartdisease'],{'chol':100})
print(query['heartdisease'])